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Remarks 

In the Office Action dated May 17, 2005, claims 49-66, in the above- 
identified U.S. patent application were rejected. Reconsideration of the 
rejections is respectfully requested in view of the above amendments and the 
following remarks. Claims 49-66 remain in this application and claims 1-48 have 
been canceled. 

Claim 49 was rejected under 35 USC §1 03(a) as obvious over Celeste, 
Ben-Bassat, and Hirel in view of Georgiou further in view of Thompson and 
Tonouchi. Applicants respectfully point out that the Pro-MP52 according to the 
present claims cannot be achieved by the process suggested in the office action. 
A pure composition according to the present invention is possible because 
proteins with Ala or Met-Ala at the N-terminus are not expressed. It is known that 
an incomplete cleavage of methionine is achieved when treating with MAP under 
both in vitro and in vivo conditions. Thus, if aminopeptidases are used after 
expression the cleavage would be incomplete even when an excess amount of 
aminopeptidases is used and thus the preparation would not have the desired 
homogeneity and the required purity. Even with the help of aminopeptidases, a 
Pro-MP52 can never be obtained which is in fact free of proteins having an Ala or 
Met-Ala at the N-terminus. In addition, in such a preparation, aminopeptidases 
would be present along with the protein. Aminopeptidases are impurities which 
would be undesirable in a medicament. A complete separation of 
aminopeptidases, which had previously been added, would not be possible. 
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Thus, the isolated protein of the present invention would not be obtained by the 
method suggested in the office action. 

As discussed above, even after treatment with aminopeptidases there 
would still be proteins with an Ala or a Met-Ala at the N-terminus present. The 
presence of non-native N-terminal methionine in recombinant pharmaceutical 
proteins can lead to problems. Applicants point out that a difference in one 
methionine residue at the N-terminal end is known to change the mobility of 
proteins in RP-HPLC and SDS-PAGE. Differences have been reported for 
recombinant ricin A chain and interleukin-2. This type of heterogeneity among 
the protein molecules can cause difficulties in the purification and 
characterization of these proteins and heterogeneity decreases the yield of pure 
products. In addition, an extra methionine at the N-terminal end of the protein 
can affect the biological activity and function of the protein, can alter its stability, 
and can make it immunogenic when administered as a therapeutic product. The 
recent debate on whether methionine-free human growth hormone (produced by 
Eli Lilly and Company) is superior to methionylated human growth hormone 
(produced by Genentech) exemplifies the importance of this issue. Applicants 
also point out that the methionine at the N-terminal end is labile to oxidation 
during the purification process, especially when treated at high pH and under 
aerobic conditions. Oxidation of the methionine residue to sulfonyl methionine 
can increase the immunogenicity of the protein. 

In addition to the above discussed disadvantages, when using 
aminopeptidases, it can be assumed that the aminopeptidase preparations 
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themselves had been isolated and purified and thus there is a risk that they are 
contaminated with further proteases which have similar purification properties. 
These contaminated proteases could attack the MP52 protein under favorable 
reaction conditions (30°C and/or 37°C overnight for the MAP). 

When MP52 proteins are present as Met-Ala-Pro-MP52, under realistic 
conditions in a manufacturing process, when using reaction conditions for 
aminopeptidases, a large part of the N-termini of Met-Ala-Pro-MP52 would not 
be accessible to the MAP due to precipitations/aggregations. An in vitro 
treatment of proteins and thus also of Met-Ala-Pro-MP52 with aminopeptidases 
requires the use of pH values around the physiologic range, i.e. around pH 7. 
The maintenance of particular pH ranges is important for enzymatic reactions so 
that an enzyme can work efficiently. The pH profile for the MAP from E.coli can 
be gathered from Figure 5 of the article Klein et al. (2003). The requirements 
necessary for the aminopeptidase reaction (pH value about 7.5) are very 
problematic for the MP52 proteins. MP52, in contrast to the aminopeptidase, is a 
highly hydrophobic enzyme, which, accordingly, has a very poor solubility under 
physiologic conditions around a neutral pH value in an aqueous solution. 
Applicants point out WO 03/043673 which deals with the homogeneous loading 
of carrier materials with 13MPs, especially human GDF-5 (MP52) in order to 
obtain devices with osteoinductive properties. Also for the loading of carrier 
materials, it is important that the protein remains in the solution as long as 
possible (cf. e.g. page 5, last sentence in the last-but-one paragraph: "An 
important precondition for this process is a sufficient solubility of the proteins in 
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the coating solution." or on page 6, bottom: "(a) providing a solution comprising 
dissolved osteoinductive protein and a buffer which keeps said protein 
dissolved... ." ). As described in this application, it is very important to keep a 
stable, suitably low pH value. A carrier material such as CaP effects an increase 
of the pH value which must be avoided by a buffer since otherwise 
aggregations/precipitations arise (cf. page 7, middle until page 8, middle). The 
pH-dependent solubility of human GDF-5 (MP52) is shown in Example 3 and 
Figures 1 and 2. It is shown that GDF-5 is only soluble in an acid range (pH 5) 
or in the alkalic range (pH > 10). In the pH range between 6.0 and 9.5 (this pH 
range is, as explained above , necessary for the aminopeptidase reaction) only 
an extremely low solubility is shown (< 5pg/ml) (cf. page 22, last two sentences 
under Example 3). 

MP52 has an extremely poor solubility under conditions wherein the 
enzymatic reaction of aminopeptidases takes place. Under conditions where 
MP52 is soluble, the MAP cannot work or can work only very inefficiently. Thus, 
in a manufacturing process and at the conditions for an aminopeptidase 
treatment, micro-precipitations and aggregations of MP52 will always occur. 
With such an aggregation/precipitation it must be assumed that a part of the N- 
termini of Met-Ala-Pro-MP52 will not be accessible to the MAP due to physical 
reasons, completely independent from the amount or the excess of MAP used. 
With regard to these N-termini it is reasonable to assume that a cleavage of Met 
would not be affected. In any case, one skilled in the art would not subject MP52 
to conditions suitable for MAP since it is known that precipitations and micro- 
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precipitations, as they occur due to poor solubility, lead to a partial or complete 
inactivation (cf. page 6, lines 4-5: "Moreover, it is to be understood that the 
osteoinductive proteins are not aggregated and partially or entirely inactivated 
due to precipitation or micro-precipitation..."). The 3-dimensional folding of the 
proteins changes due to precipitation. 

Tonouchi does not make any statements regarding the pH-value but refers 
to an article of Yaron (1971). Yaron indicates that the optimal pH value is at 8.6 
(last sentence, page 529). Figure 1 A on page 530 shows the pH profile of 
aminopeptidase P. At the end of the first paragraph on this page, Yaron states 
that "This maximum was observed at all pH values at which activity could be 
tested (7.0-9.2)". The pH range of 7.0-9.2, wherein the activity of 
aminopeptidase P can be tested, is always in a range where MP52 has a poor 
solubility (sufficient solubility of MP52 occurs only above pH 10) and therefore a 
part of the N-termini of Met-Ala-Pro-MP52 would not be accessible to the MAP. 

In view of the above discussion, it is clear that during the treatment of Met- 
Ala-Pro-MP52 with aminopeptidases, some of the methionine would not be 
cleaved even when using an excess of MAP. Thus, the Pro-MP52 preparation 
obtained using the process suggested in the office action would be contaminated 
with an excess of MAP and aminopeptidase P, with Met-Ala-Pro-MP52 and 
possibly also with Ala-Pro-MP52. Thus, such a preparation clearly differs from 
the pure Pro-MP52 preparation according to the present claims. 

Applicants also contend that it would not have been obvious from the prior 
art to produce Pro-MP52 at all. In other words, at the time of the present 
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invention the skilled artisan would not have been motivated to produce MP52 as 
a homogeneous material since there is no suggestion in the prior art that Pro- 
MP52 could be obtained in a homogeneous form and that it would not be further 
decomposed by aminopeptidases. The use of aminopeptidases would lead to 
totally different mixed products and not to homologously isolated proteins as in 
claim 49. In view of the above discussion, applicants request that this rejection 
be withdrawn. 

Claims 49 and 50 were rejected under 35 USC §1 03(a) as obvious over 
Celeste, Ben-Bassat, and Hirel in view of Georgiou further in view of Thompson 
and Tonouchi and further in view of Hotten and Cerletti. Claims 49-60, 63 and 66 
were rejected under 35 USC §1 03(a) as obvious over Celeste, Ben-Bassat, and 
Hirel in view of Georgiou further in view of Thompson and Tonouchi and further 
in view of Hotten and Cerletti and Neidhardt. Claims 49-51 , 64, 65 and claims 
52-60, 63 and 66 were rejected under 35 USC §1 03(a) as obvious over Celeste, 
Ben-Bassat, and Hirel in view of Georgiou further in view of Thompson and 
Tonouchi and further in view of Hotten and Cerletti and Neidhardt further in view 
of Hotten A and Chen. Claims 49-51, 61 and 62 were rejected under 35 USC 
§1 03(a) as obvious over Celeste, Ben-Bassat, and Hirel in view of Georgiou 
further in view of Thompson and Tonouchi and further in view of Hotten and 
Cerletti and Neidhardt further in view of Ron and Avis. Claims 49-51 and 60 
were rejected under 35 USC §1 03(a) as obvious over Celeste, Ben-Bassat, and 
Hirel in view of Georgiou further in view of Thompson and Tonouchi and further 
in view of Hotten and Cerletti and Neidhardt further in view of Oppermann. As 
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discussed above, a homogeneous Pro-MP52 cannot be obtained when 
aminopeptidases are used because the enzymes are not 1 00% efficient. 
Therefore, applicants contend that the limitation in claim 49 that proteins 
according to SEQ ID NO:1 with either a) an ala or b) a Met-Ala at the N-terminus 
are not expressed or present distinguish the claimed invention from the cited 
prior art. In view of the above discussion applicants request that these rejections 
be withdrawn. 

Claims 49-66 were rejected under the judicially created doctrine of 
obviousness type double patenting as unpatentable over claims 1-7 of copending 
application no. 10/365,231. A terminal disclaimer is being prepared and will be 
filed shortly. 

Applicants respectfully submit that all of claims 49-66 are now in condition 
for allowance. If it is believed that the application is not in condition for 
allowance, it is respectfully requested that the undersigned attorney be contacted 
at the telephone number below. 
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In the event this paper is not considered to be timely filed, the Applicant 
respectfully petitions for an appropriate extension of time. Any fee for such an 
extension together with any additional fees that may be due with respect to this 
paper, may be charged to Counsel's Deposit Account No. 02-2135. 



Respectfully submitted, 

Monica Chin Kitts 

Attorney for Applicants 

Registration No. 36,105 

ROTHWELL, FIGG, ERNST & MANBECK, p.c. 

Suite 800, 1425 K Street, N.W. 

Washington, D.C. 20005 

Telephone: (202)783-6040 
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We have performed a computational study of different 
protomeric states of the methionine aminopeptidase ac- 
tive site using a combined quantum-mechanical/molec- 
ular mechanical simulation approach. The aim of this 
study was to clarify the native protonation state of the 
enzyme, which is needed for the development of novel 
irreversible inhibitors that can possibly be used as an- 
tiangiogenic and antibiotic drugs by virtual screening 
and other drug design methods. The results of the sim- 
ulations indicated that two protonation states are pos- 
sible without disturbing the overall geometry of the ac- 
tive site. We then verified experimentally the presence 
of the two protonation states by studying the substrate 
hydrolysis and inhibitor binding reactions at different 
pH values and come to the conclusion that one of the 
protomeric states is relevant for inhibitor binding, 
whereas the other is relevant for substrate hydrolysis. 
This result has implications for the development of 
other inhibitors of this class of enzymes and adds a new 
perspective to the pharmacological properties of the 
antiangiogenic drug fumagillin, which is an irreversible 
inhibitor of the human methionine aminopeptidase 
type EC. 



Methionine aminopeptidase s (MetAPs) 1 play a central role 
for in vivo protein synthesis as they remove the starter methi- 
onine from newly synthetized proteins. The natural product 
fumagniin, an epoxide, covalently modifies one of the active- 
site histidines in the eukaryotic methionine aminopeptidase II 
(MetAP-II) (1-4) and other MetAPs (of. Fig. 1). Fumagillin 
inhibits the growth of vessels in tumors, and a derivative of the 
compound has been evaluated in clinical trials as an anticancer 
drug. The antiangio genie effect of fumagillin and other inhibi- 
tors of MetAP-II have been attributed to the inhibition of the 
Ets-1 transcription factor expression and the activation of the 



* The costs of publication of this article were defrayed in part by the 
payment of page charges. This article must therefore be hereby marked 
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p53 pathway (5, 6). Besides from being anticancer drug targets, 
MetAPs have the potential to become the target proteins of 
antibacterial substances because the MetAP functionality is 
essential for cell growth and bacteria possess only one of the 
two known MetAP subtypes (7). 

MetAPs are metal-dependent enzymes. It is not clear which 
metal activates the MetAPs in vivo. For" in vitro experiments, 
cobalt is commonly used because it activates all known MetAPs 
and the cobalt-substituted enzymes are usually the most active 
(8). Zinc and iron(II) have also been shown to activate some 
MetAPs (9, 10). The metal-chelating residues in all known 
MetAPs are two glutamates, two aspartates, and one histidine. 
The geometric arrangement of these residues is practically 
identical in all of the MetAP x-ray structures (11, 12). Several 
three-dimensional structures of MetAPs have been determined 
by x-ray diffraction methods including the structure of human 
MetAP-II with a covalently bound fumagillin molecule (13-15). 

Density functional theory (16, 17) has long been used to 
study the reactivity of organic molecules and recently began to 
find its way into biochemistry (18). We present here an exam- 
ple in which the application of density functional theory to a 
biochemical problem led to a hypothesis, which was subse- 
quently verified by experiments. Our results do not only ex- 
plain a couple of biochemical phenomena related to MetAPs but 
also demonstrate the impact that modern quantum-chemical 
methods can have on the study of biological systems where they 
can help fill the gap between mere conjectures and physical 
reality. 

Our motivation for examining the MetAPs from a theoretical 
point of view stems from the highly selective, irreversible in- 
hibition mechanism of fumagillin and related epoxides. Many 
enzyme inhibitors that form a covalent adduct with their target 
proteins (e.g. the /3-lactam antibiotics and acetyl salicylic acid) 
are important drugs. We believe that the nimagillin/MetAP 
example is a good test case for theoretical methods that aim at 
rationalizing the development of covalent enzyme inhibitors, 
because a large amount of high quality x-ray structural data, 
also with bound inhibitors, has been published over the last few 
years. 

As a first step toward an understanding of the catalytic and 
inhibitor-binding mechanisms of MetAPs, we studied different 
protonation states of the active site water molecules by means 
of molecular dynamics simulations in a combined quantum- 
mechanical/molecular-mechanical (QM/MM) framework. In the 
1.9 -A resolution, x-ray structure of Escherichia coli MetAP 
(Protein Data Bank code 2mat), one water molecule (or hydrox- 
ide ion) is bridging the two cobalt ions and another water 
molecule is bound to the cobalt ion that is not coordinated to the 
histidine. Three different protomeric states were examined (see 



47862 



This paper is available on line at http://www.jbc.org 



Protonation States of Methionine Aminopeptidase 



47863 



Fig. 2) as follows: one with two bound water molecules, one 
with a bridging hydroxide ion and a water molecule, and one 
with two hydroxide ions. For the quantum -mechanical part of 
the system, the ab initio molecular dynamics method described 
by Car and Parrinello (CPMD) (19) was used. Because an 
accurate treatment of the open-shell cobalt(II) ion is very de- 
manding in the framework of CPMD, we decided to use zinc(II) 
ions as the active-site metals in the CPMD simulations. The 
fundamental problem with cobalt is the uncertainty of its spin 
state and the multiplicity of the bi-cobalt system. A computa- 
tional study of a di-cobalt system would therefore require the 
simulation of eight different electronic configurations that is 
not feasible with CPMD and currently available supercom- 
puter resources. Given the fact that at least some MetAPs have 
been shown to be active with zinc and considering that a 
thorough theoretical analysis of the zinc- and cobalt-substi- 
tuted truncated MetAP active sites did not reveal major geo- 
metrical differences between a di-zinc and a di-cobalt system 
(20), we assume that the replacement of cobalt by zinc does not 
affect the validity of our results. Furthermore, preliminary 
CPMD simulations on a di-cobalt MetAP active site with one 
water and one hydroxide ion did not show geometric or dynamic 
differences to a di-zinc system (see Supplemental data). 

EXPERIMENTAL PROCEDURES 

Theoretical Methods — The density functional theory-based Car -Par- 
rinello (19) molecular dynamics program CPMD, version 3.5, was used 
for the QM/MM simulations of MetAP (21-24). AMBER94 force-field 
parameters were used for the MM part (25, 26). All of the systems were 
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Fig. 1. Binding of fumagillin in the MetAP active site. Me, metal 
ion. The numbering is for E. coli MetAP. 



fully hydrated with TIP3P (three-point transferable intermolecular po- 
tential) water molecules in a periodic box and equilibrated by extensive 
classical-mechanical molecular dynamics (MD) simulations with AM- 
BER94 prior to the ab initio MD simulations. Visualization was done 
with VMD (27). The quantum parts included the five metal-binding 
amino acids (truncated to acetate ions and imidazole, respectively; free 
valences were capped by adding hydrogen atoms), the two zinc ions, and 
the two water/hydroxide molecules coordinating to them. The QM parts 
were minimized (annealed to a temperature below 0.1 K) before the 
start of the production MD phase. CPMD parameters were as follows: 
isolated system calculations; gradient-corrected exchange-correlation 
functionals due to Becke (28) and Lee et at. (29); plane waves basis set; 
kinetic energy cutoff: 70 Ry; soft norm-conserving Troullier-Martins 
(30) pseudopotentials for the core electrons; time step (6 atomic units); 
Nose-Hoover thermostat at 300 K, coupling frequency of 500 cm" 1 ; and 
fictitious electron mass (800 a.u.) (28-32). The size of the orthorhombic 
QM box was 13.2 X 15.9 X 17.2 A, corresponding to a minimum image 
distance of 6.4 A. Simulations with cobalt were performed with a kinetic 
energy cut-off of 90 Ry and Becke (28) and Perdew (33) exchange- 
correlation functionals (33, 34). 

Experimental Methods — E. coli MetAP was purified from an overpro- 
ducer strain that was kindly provided by Drs. Lowther and Matthews. 
Assays were performed in 96-well microtiter plates using the method of 
Yang et al. (35) with Met-Gly-Met-Met as substrate and an enzyme 
concentration of 12 nM. The fumagillin concentration in the inhibition 
experiments was 66 /am, which is slightly below the IC 50 of fumagillin at 
pH 7 under otherwise identical conditions. Tris/maleate buffer was used 
for the incubation of MetAP with and without fumagillin at pH 5—8.5 
(volume of 30 /d, 15 min, 37 °C). For the ensuing Met-Gly-Met-Met 
hydrolysis and detection reactions, the pH was adjusted to pH 7.5 by 
adding 170 fil of 4-fold concentrated Tris/maleate buffer. The fluores- 
cent reaction product (resorufin) was measured using a Wallac micro- 
titer plate reader. All of the experiments were performed in triplicate. 
Detailed descriptions of the experiments are given in Supplemental 
material. 

RESULTS 

Simulations — In the CPMD simulations, the system with the 
two water molecules coordinated to the zinc ions showed a 
pronounced movement of the coordinating water molecules 
away from the x-ray structure of E. coli MetAP that has one 
water molecule located between the metals. The system 
evolved to a structure in which there was no more "bridge" 
between the metals and one water molecule started to form a 
strong hydrogen bond to a carboxylic acid. The simulation was 
stopped after —5000 steps of simulation (0.85 ps), because the 
geometry remained essentially stable. In contrast, the struc- 
ture with one bridging hydroxide was stable during the whole 
simulation time (~10,000 steps), showing only the expected 
thermal motions and a twist of one carboxylic acid group. Fig. 
3 gives a comparison of the active site geometries of the bi- 
water and the water-hydroxide system with the crystal struc- 
ture. Hydrogens have been omitted for clarity. The final state of 
the bi-hydroxide system is not shown, because it deviates 
strongly from the crystal structure and a comparison with the 
other systems would not make any sense. After only —1000 
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Fig. 2. A-C, the three different protomeric states for which CPMD simulations were performed. 
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Fig. 3. Final structures of the 
CPMD simulations of the bi-water 
and the water-hydroxide systems in 
comparison with the x-ray structure. 

Hydrogens and the surrounding residues 
of the enzyme are not shown for clarity. 
The x-ray-structure is shown in CPK 
(Corey-Pauling-Koltun) colors, the bi-wa- 
ter system in yellow and the water-hy- 
droxide system is in green. Me, metal ion; 
Wat, water. 





Fig. 4. Final structure of the active site with one hydroxide ion and one water molecule coordinated to the metals (Mel, Me2, 
purple spheres). The HOMO-1 orbital located at the bridging hydroxide ion is also shown (red and blue clouds). The cut-off for the visualization 
of the HOMO-1 electron density was 0.1 e/au 3 (electron per cubic atomic units). 



steps of simulation, large deviations from the crystal structure 
were observed in the bi-hydroxide system. After 3000 steps, the 
metal-coordinating residues were considerably displaced from 
the crystal structure positions and the simulation was stopped. 

It becomes clear from Fig. 3 that both the bi-water and the 
water-hydroxide system are intrinsically stable and quite sim- 
ilar to the 2 mat x-ray structure. The root mean square devia- 
tions of the active site heavy atoms in comparison with the 
x-ray structure (hydrogens were not used because they are not 
present in the x-ray structure) are 0.81 A for the bi-water and 
0.85 for the water-hydroxide system. The main changes are as 
follows. 1) The bridging water leaves the position between the 
two metals in the bi-water simulation. 2) In the water-hydrox- 
ide simulation, the bridging hydroxide remains in place and 
one of the metal-coordinating carboxylic acid groups twists to a 
position that is perpendicular to the x-ray structure. Neverthe- 
less, this residue remains a binding partner for one of the metal 
ions. 

A picture of the final state of the simulation is shown in Fig. 
4. This figure also shows the HOMO-1 orbital, which is local- 



ized at the hydroxide ion, the species that very likely acts as the 
nucleophilic agent in the substrate hydrolysis reaction. The 
HOMO orbital, which is —0.5 eV higher in energy, is localized 
at one of the carboxylic acid groups. 

Experiments — The pH profiles of the inhibitor binding and 
substrate hydrolysis reaction are given in Fig. 5. It is evident 
that the pH optimum of the substrate hydrolysis reaction is at 
pH 7.5, whereas the binding of fumagillin becomes increasingly 
favored at lower pH. 

DISCUSSION 

The results presented above allow the following conclusions. 
1) In the x-ray structure of E. coli MetAP (2mat), the bridging 
water molecule is most probably deprotonated. However, an 
active site with two (fully protonated) water molecules coordi- 
nated to the metal ions is also stable, albeit with a different 
coordination geometry. The x-ray structure of the E. coli 
MetAP (2mat) was determined at a pH larger than 7, whereas 
several other MetAP structures were determined at a more 
acidic pH. (human MetAP-H, see Ref. 15), The pH in the crys- 
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Fig. 6. Mechanisms of the fumagil- 
lin binding and substrate hydrolysis 
reactions in MetAP. Only the initial 
steps of the substrate hydrolysis are 
shown (see the literature references in 
the text for more details regarding the 
MetAP substrate hydrolysis mechanism). 



0 etc. 



tallization medium can not be deterinined exactly, but the 
experimental (buffer) conditions described by Liu et al. (15) 
indicate that the pH was below 7. Examining one of these x-ray 
structures (lbn5), we found that no bridging water molecule is 
present between the two metals. Taken together, the x-ray 
structure and the results of the CPMD simulations indicate 
that two protonation states are possible without disrupting the 
active site geometry. At a more basic pH, one water molecule is 
deprotonated and bridging the two metal ions. From the elec- 
tronic structure calculations, this protonation state is expected 
to be relevant for the catalytic process. Two water molecules 
are present in the active site at more acidic pH, each one 
coordinating a different metal ion. 

The covalent binding of fumagillin to MetAP requires proto- 
nation of the epoxide oxygen. Proton donation from an active 
site group, in particular from one of the active site water 



molecules, is more Ukely in the protonation state with two 
water molecules (i.e. at more acidic pH). 

In other words, our calculations and the analysis of the 
crystal structure indicate that the active site of binuclear 
E. coli MetAP is characterized by the presence of a water - 
molecule with a piT c of ~7. This is at the lower border of the 
pK a values reported for waters coordinating to zinc ions in 
organic complexes (36). The binding of fumagillin is expected to 
be favored at more acidic pH, whereas a more basic pH would 
be required for catalysis. This finding opened the way to an 
experimental validation of the theory of which we determined 
the pH-profile of the fumagillin inhibition reaction and the 
enzymatic activity of MetAP (see Fig. 5). It is obvious that the 
binding of the inhibitor is favored under more acidic conditions 
as compared with those that are optimal for the substrate 
hydrolysis reaction. This observation provides strong substan- 
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tiation to our interpretation of the QM/MM CPMD simulation 
results. 

In summary, the mechanisms of the fumagilUn binding and 
the initial step of the substrate hydrolysis reactions are shown 
in Fig. 6. The substrate hydrolysis reaction starts with the 
attack of the metal-bound hydroxide ion to the carbonyl carbon 
of the scissile peptide bond. This mechanism differs slightly 
from the one that is given in the recent literature (12) where 
the deprotonation of the metal-bridging water occurs after the 
substrate has been bound in the active site. 

Recent theoretical work on the reaction mechanism of imi- 
penem binding to the di-nuclear zinc-/3-lactamase (37) and the 
substrate hydrolysis mechanism of the di-nuclear bovine lens 
leucine aminopeptidase (38) indicate that the nucleophilic at- 
tack on the substrate originates from a "terminally* bound and 
not a bridging hydroxide ion. The nucleophilicity of a hydroxide 
ion that is complexed to one metal ion is higher than that of a 
bridging hydroxide. Considering these results, we suggest a 
mechanism (Fig. 6, bottom) in which a bridging hydroxide (as in 
the crystal structure) is found in the resting state of the enzyme 
and the binding of the substrate causes a rearrangement of the 
complexation pattern. Subsequently, the amide carbonyl car- 
bon is attacked by the terminally bound hydroxide ion, which 
has an increased nucleophilicity as compared to the bridging 
hydroxide ion in the resting state. 

It may be argued that the acidic pH leads to the protonation 
of an active site residue other than the metal-coordinating 
water and that fumagillin binding is therefore facilitated by a 
modification of the non-covalent complex that precedes the 
binding reaction. However, the only groups that could readily 
be protonated within the studied pH range would be His-79 and 
His-178. His-79 is the binding partner of fumagillin. A proto- 
nated and therefore positively charged His-79 would certainly 
be much less prone to act as a nucleophile on fumagillin, and it 
can thus be assumed that His-79 remains neutral. Therefore, if 
one assumes that His-79 is not being protonated, it is also very 
unlikely that His- 178 is protonated because the chemical sur- 
rounding (which could influence the pK a of the histidine imid- 
azole) is quite similar for both histidines. Furthermore, litera- 
ture values for the histidine side chain pK a range from 6 to 6.4, 
which is below the value of ~7 for the protonation event in the 
MetAP active site. 

The results of the CPMD simulations presented here do not 
provide quantitative free energy data of the different proto- 
meric states. It would certainly be very useful to obtain such 
data, which could be used to calculate the pK a of metal-bound 
waters and other ionizable groups in the active site. However, 
one must be aware of the fact that the prediction of pK a values 
would require very exact free energy data of different proto- 
meric states and a very complete sampling of conformational 
space. Even in the classical-mechanical simulation framework 
(where the computational cost of simulations is smaller by 
several orders of magnitude and a much more complete explo- 
ration of conformational space is possible), free energies are 
notoriously difficult to predict (39). Thus, we think that a 
meaningful prediction of free energy differences and pK a values 
by ab initio simulation methods is currently not feasible. The 
situation may be somewhat different when less complex sys- 
tems such as pure water or small solvated ions are considered. 
In such cases, present computer capabilities can be sufficient to 
adequately determine pK a values by ab initio MD simulations. 

Recent experimental work (40) has shown that E. coli MetAP 
may under certain circumstances function as a mononuclear 
enzyme, i.e. with only one metal ion coordinating to the active 
site residues. The metal content and the identity of the metal in 
the MetAP enzymes in vivo are not known. With the sole 




Fig. 7. De-pro tonation of a 0-hydroxynitrile ester, subsequent 
ester hydrolysis, and formation of acrylonitrile. Acrylonitrile then 
binds to micleophiles in the active site, e.g. histidine or cysteine 
residues. . 

exception of the publication mentioned above, all of the exper- 
imental data of E. coli and yeast and human MetAPs, including 
all of the published x-ray structures, have been obtained with 
bi-nuclear enzymes because enzymes with lower metal content 
had only marginal or no catalytic activity under the respective 
experimental conditions. Therefore, we are confident that it is 
reasonable to assume the presence two metal ions in the active 
site as we did in our simulation system. 

Fumagillin and its congeners inhibit the growth of vessels in 
tumors. Considering the fact that the extracellular pH in tu- 
mors is more acidic than in normal tissues, (41) and assuming 
that the endothelial cells of vessels in tumors are under the 
influence of the acidic extracellular tumor pH, one may reason 
that the pH profile of the fumagillin-MetAP binding reaction 
leads to selectivity (or targeting) of the fumagillin effect to 
tumor vessels. 

The results presented here have several implications for the 
rational design of MetAP inhibitors. First, virtual screening 
(docking/scoring) procedures should preferably be made on 
MetAP active sites with a net charge of - 1 if inhibition at 
neutral pH is desired (e.g. for MetAP inhibitors with antibac- 
terial activity). The active site net charge has a pronounced 
influence on the binding energetics of ligands in virtual screen- 
ing and in physical reality, because electrostatic interactions 
are the strongest driving force for intermolecular attraction 
and repulsion. Second, irreversible inhibitors different from 
fumagillin should preferably be compounds that are activated 
by a nucleophilic attack (coming from the hydroxide ion) or by 
deprotonation. Such inhibitors would be much more active at 
physiological pH than fumagillin. For instance, the deprotona- 
tion of /3-hydroxynitrile esters (cf. Fig. 7) causes cleavage of the 
ester bond and subsequent formation of acrylonitrile, which 
can then bind (as a ''Michaelis-like* electrophile) to one of the 
cysteine residues near the MetAP active site. The sequence of 
reactions would be based on a "basic activation" and not acidic 
activation as with fumagillin. 

In conclusion, we have clarified the protonation behavior of 
the MetAP active site and made an experimental proof of a 
hypothesis that was generated on the basis of quantum- 
mechanical calculations. 

Acknowledgments — We thank Prof. R. Hart maim at Saarland Uni- 
versity for continued support for this project and the Deutsche Forsch- 
ungsgemeinschaft for financial support to C. Klein (KL 1356). CSCS in 
Manno, Switzerland and HLRS in Stuttgart, Germany kindly provided 
CPU time for this project. The group of Prof. R. Bernhard at the 
University of the Saarland helped with the expression of E. coli MetAP 
from a strain that was kindly provided by Profs. B. Matthews and 
W. T. Lowther. We are also grateful to one of the reviewers, who 
suggested us to consider a terminally bound hydroxide ion as the 
nucleophilic species in the substrate hydrolysis reaction. 

REFERENCES 

1. Ingber, D., Fujita, T. f Kishimoto, S., Sudo, K., Kanamaru, T., Brem, H., and 

Folkman, J. (1990) Nature 348, 555-557 

2. Griffith, E. C, Su, Z., Turk, B. E., Chen, S., Chang, Y. H., Wu, Z., Biemann, K., 

and Liu, J. O. (1997) Chem. Biol 4, 461-471 

3. Sin, N., Meng, L. t Wang, M. Q. W., Wen, J. J., Boramann, W. G., and Crews, 

C. M. (1997) Proc. Natl. Acad. Sci. U. S. A. 94, 6099-6103 

4. Griffith, E. C, Su, Z., Niwayama, S., Ramsay, C. A., Chang, Y. H., and Liu, 

J. O. (1998) Proc. Natl. Acad. Sci. U. S. A. 95, 15183-15188 



Protonation States of Methionine Aminopeptidase 



47867 



5. Wernert, N., Stanjek, A., KiriaMdis, S., Hiigel, A., Jha, H. C, Mazitschek, R., 

and Giannis, A. (1999) Angew. Chem. Int. Ed. Engl. 38, 3228-3231 

6. Zhang, Y., Griffith, E. C, Sage, J., Jacks, T., and Liu, J. O. (2000) Proc. Natl. 

Acad. Sci. U. S. A. 97, 6427-6432 

7. Li, X., and Chang, Y. H. (1995) Proc. Natl. Acad. Sci. U. S. A, 92, 12357-12361 

8. Ben-Bassat, A., Bauer, K, Chang, S. Y., Myambo, K., Boosman, A, and Chang, 

S. (1987) J. Bacteriol. 169, 751-757 

9. D'Souza, V. M., and Holz, R. C. (1999) Biochemistry 38, 11079-11085 

10. Walker, K. W., and Bradshaw, R. A (1998) Protein Sci. 7, 2684-2687 

11. Lowther, W. T., and Matthews, B. W. (2000) Biochim. Biophys. Acta 1477, 

157-167 

12. Lowther, W. T. t and Matthews, B. W. (2002) Chem. Rev. 102, 4581-4608 

13. Roderick, S. L., and Matthews, B. W. (1993) Biochemistry 32, 3907-3912 

14. Lowther, W. T., Orville, A. M., Madden, D. T., Lim, S., Rich, D. H., and 

Matthews, B. W. (1999) Biochemistry 38, 7678-7688 

15. Liu, S., Widom, J., Kemp, C. W., Crews, C. M., and Clardy, J. (1998) Science 

282 1324—1327 

16. Kohn, W., and Sham, J. (1965) Phys. Rev. A 140, 1133-1138 

17. Hohenberg, P., and Kohn, W. (1964) Phys. Rev. B 136, 864-871 

18. Carloni, P., Rothlisberger, U., and Parrinello, M. (2002) Acc. Chem. Res. 35, 

455-464 

19. Car, R., and Parrinello, M. (1985) Phys. Rev. Lett. 55, 2471-2474 

20. Jorgensen, A T., Norrby, P.-O., and Liljefors, T. (2002) J. Comput. Aided Mol. 

Des. 16, 167-179 

21. Hutter, J., Alavi, A., Deutsch, T., Bernasconi, M., Goedecker, S., Marx, D., 

Tuckerman, M., and Parrinello, M. (1995-2003) CPMD, Max-Planck Insti- 
tute fur Festkbperforschung, Stuttgart, Germany and IBM Zurich Research 
Laboratory, Zurich, Switzerland 

22. Laio, A., VandeVondele, J., and Rothlisberger, U. (2002) J. Chem. Phys. 116, 

6941-6947 

23. Laio, A., VandeVondele, J., and Rothlisberger, U. (2002) J. Phys. Chem. B 106, 

7300-7307 

24. Colombo, M. C, Guidoni, L., Laio, A, Magistrate, A, Maurer, P., Piana, S., 

Rohrig, U., Spiegel, K, Sulpizi, M., VandeVondele, J., Zumstein, M., and 



Rothlisberger, U. (2002) Chimin 56, 11-17 

25. Cornell, W. D., Cieplak, P., Bayly, C. I., Gould, I. R., Merz, K. M., Ferguson, 

D. M., Spellmeyer, D. C, Fox, T., Caldwell, J. W. ( and Kollman, P. A (1995) 
J. Am. Chem. Soc. 117, 5179-5197 

26. Case, D. A, Pearlman, D. A, Caldwell, J. W., Cheatham III, T. E., Ross, W. S., 

Simmerling, C. L., Darden, T. A, Merz, K M., Stanton, R. V., Cheng, A L., 
Vincent, J. J., Crowley, M., Tsui, V., Radmer, R., Duan, Y., Pitera, J., 
Massova, I., SeibeL G. L., Singh, U. C, Weiner, P., and Kollman, P. A 
(1999) AMBER6, University of California, San Francisco, CA 

27. Humphrey, W., Dalke, A, and Schulten, K (1996) J. Mol. Graph. 14, 33-38 

28. Becke, A. D. (1988) Phys. Rev. A. 38, 3098-3100 

29. Lee, C. L., Yang, W., and Parr, R. G. (1988) Phys. Rev. B. 37, 785-789 

30. Troullier, N., and Martins, J. L. (1991) Phys. Rev. B. 43, 1993-2006 

31. Nose, S. (1984) Mol. Phys. 52, 255-268 

32. Barnett, R. N., and Landman, U. (1993) Phys. Rev. B. 48, 2081-2097 

33. Perdew, J. P. (1986) Phys. Rev. B. 33, 8822-8824 

34. Rovira, C, Kunc, K, Hutter, J., and Parrinello, M. (2001) Inorg. Chem. 40, 

11-17 

35. Yang, G., Kirkpatrick, R. B., Ho, T., Zhang, G. F., Liang, P. H., Johanson, 

K O., Casper, D. J., Doyle, M. L., Marino, J. P., Thompson, S. K., Chen, 
W. F., Tew, D. G. f and Meek, T. D. (2001) Biochemistry 40, 10645-10654 

36. Bertini, I., and Luchinat, C. (1994) in Bioinorganic Chemistry (Bertini, I., 

Gray, H. B., Lippard, S. J., and Valentine, J. S., eds), pp. 37-106, University 
Science Books, Sausalito, CA 

37. Suarez, D., Diaz, N., and Merz, K M., Jr. (2002) J. Comput. Chem. 23, 

1587-1600 

38. Schiirer, G., Horn, A. H. C, Gedeck, P., and Clark, T. (2002) J. Phys. Chem. B. 

106, 8815-8830 

39. Reinhardt, W. P., Miller, M. A, and Amon, L. M. (2001) Acc. Chem. Res. 34, 

607-614 

40. Cosper, N. J., D'souza, V. M., Scott, R. A, and Holz, R. C. (2001) Biochemistry 

40, 13302-13309 

41. Gillies, R. J., Raghunand, N., Karczmar, G. S., and Bhujwalla, Z. M. (2002) J. 

Magn. Reson. Imaging 16, 430-450 



■si ^ ,j; o 
= o ~ — 



CD 



-1 ^ 3} S 



o 



'5 « ^ = 



o £ cu ^ 53 



r ^ ic 

a; jd_ Jr; 
a; > "~ " 



S I -i g e e 

*~ ! — — ! 



o 



S -5 °o I * 

I S o s 

m 0 " -a 

« o £ ,-3 



- 1 - J <u _ F p 
— CD 

- r~ o ^ 

^ ^ *3 

a: H o ^ 

- o cd 0 ■ 

e/2 D; oj 

^ ~ g. - 



C- 

£ O O 

C3 r "r " M 

*^+- o> 

i * 3 
.s " 

cc ^ ~ 
a: P — ' 
o ~ 

c o c: 

:Ss 

£ £ -2 



^ * o 



2 « 



3= C3 

2 rt " 



-3 « 0 % -§ 



"3 O 

P. 



o CO *H \ 

CM O G _2 

CD 07 ' 3^ — 

^ — ^ ^ 



> 3 o d5 ? 



c3 S co 10 -c**: 
CD -^J O CD 

•c: '3 t=i cd 



cu - ;r; 



.E 0 — 

w ° p 0 ^ 2 

£ tn ~ CD ^ 

« c3 Q (U CM rO 

Oh - ^ 

zz CO 

(73 CO 



o _ o 

'•5 Of) 5 

c ? ^ 

^ _2 ^ 



t-. o 
<u o 



^ S 3 ^ 



■-; O 



CO 

a 



o 

V7 

33 . 

CD ^ 



x 1 a; c -7J 



o a: o 



o 



c — 



•5 0 S 5 - ^ 

7J ^ ^ Si Ci ^ 

— d 5 5? o & 

a i ^ s:^ 



. — ; cD 
o « ^ 



^ O 

~ C 

— 1 

0 5 3 55 



"5.ti S 

C3 ^ _ 



b « o 

<x *— o 

0 ^ ^ 
^ ^ cd 

,3 g 



O o 
CD 

-| 1 

3 ^ 
— . o 



CD 



X w 



i 5 



^ v -w 



f <Z> CZ z-, 



° 11 { !i 



P =3 ^ 



~- ^ £: >- 



10 
oc- 



53 o 
S a- 



2 — .2 

O - — r/; C_ c: 

^ ^ ^ £ C3 



V5 <r o ^ 

' t! ^ 

a) F- a> — ■ 

rj" rt ■ a_- LIZ r: 



« ^ b: « 









a) 


0 










0 







a 
O 



-t3 ^ X 
§ -2 

O ^ CD ^ 



"3 ^ CD 



« s a Ch 



I— i . 

o — ■ 
O = 



o T? 

^ X 



2 5 i .5 -? 



o ,S 



c3 "m 
O 



= s ~ g. 

S? <- 2 ^~ 

■3 § o 2 

" « § E- 

— ; cj 0 

■4.3 (j^J p> 

^ id 2 ^ c 

O r/2 C 



S ^ M ~ 



CO 



0 --a 



— =3 

2 a 

CD * 

^3 =5 C 



j- CO 

,iZL 

CD O 



"■5 " H o 



.2 



cs O (L 'r c 



0 g"^ 3 

™ & * CD 

CO DC- > — * 

v- c5 ^ U< 



O 

> 

O o 
& O 
o ^1 



CD CD cd 

■4-3 Z~l 

ci O O 
C! 



o 

C 



2 



2 



•J2 
fx! 



s g c i ^ 



a) ^ ^ 



,0 £| 33 cd es ~ 

.= "H I g I g -a 

— 1=3 ^ ^ - ^""S 

3 £ ~ S > — _ 0 

- - ^ = £ § 2 
, c ^ 2 



r- o ^3 — ~ 



a 0 



~ to a> f — 
'■^ c ^ i 

^ J 5 1 ^ 



p 
PU 
cr 

CD 
CJ 

a. 



S i § 3 

cd -r ^ 

j-. <U o 

o « ^ 7^ 

o C 4O 5 



~ S 2 
-• = d, "g 



r3 «; 

CD ^ 

a> ^-t . — ; 



^1 ^2 



CD 1 •— 
■xi co -r" IT 

0 a> o s 



t=JD 



S O o 

0 ^ 



C3 



O CO 

*- w cd ^n; c 

K .2 ^ M ^ 

05 - as 



co 

CD ^ 



.2 g ^ 

us § 



-^f CD 

CO 



c 
a 



5 'TS 



AO o 1 — ' 
^3 o 



.5 ? O 
O 



c: ^ W ^ 
^ 2 ^ w| 



!= f*=i O 
a; ' — ' CD 

O.' CO 



^ ^ ^ > I 3 x - o 



[S3 



p 2 ^ 

co 22 0 



^ c^i ^ P 

QJ 



S - P 



OJ Cv OJ 
— ° ^ 



CD 



5 ^ o ^ 



S3 CD- 



-c: '-c: c? 
« .2 



? © "3 

C o cj 

O CD > 

f^O CO , 

- - 



O 



^- c3 O 



H ™ ^ O 

^ CD 

-c; cd 



■3 g J §»- 

-7J ^ O 



C3 



T- -3 O 



QJ ^ ^ 



CD , 

» !fl C ^ 



— ' <~n — _ 



c dV *— 



CD 
CD 



<D PC 



— 7. hp _ 

~~ § iH p E 



ro 



.2 £ 



cu 


O 








O 


w net 1 


cell 


C3 


' VJ 


yr: 
*o 


CD 




0 

CD 

LCD 


a> 
ro 
ct 



CD « — _ ? 



3: 10 

~ CM 



* I 



p — 



I 



Procedure 



Cell-free extract. 

Heat precipi lation 

Ammonium sulfate fractionation 

Acetone fractionation 

Ca-phosphate gel 

DEAE-ccllulose 

Gel filtration 



TABLE I 

PI7IUFICATIOK OF AAnNOI»EFTIDASE-P c 



Volume Activity 6 Protein 6 Specific activity 

(ml) (units/ ml) Total' unite fmg/ml) (units/mg) 



Yield 



6,300 
6,300 
1,000 
1,000 
2,000 
29 
107 



.3 
.23 



10. 
0. 
4. 
197 
42. 



14,500 
•14,100 
10,500 
9 ..500 
8,070 
5,720 
4,560 



16.9 
11.3 
22 



6 
78 
23 
62 
37 



0.136 
0 . IBS 
0.465 
2.52 

17,3 

42.5 
116 



100 
96.7 
72 
05 
56 

39.3 
31.4 



Purification 
factor 



1 

1.45 
3.4 
IS. 4 

127 

313 

S55 



B Figures are given for 1300 g frozen wet cells of E. coli B. 

» The activity was determined as described in the text. One unit of activity is defined as the amount of enzyme ^hich produces 3 mg 

proline per milliliter of incubation solution under the conditions of the above-desenbed assay. 
* Protein concentration was determined by the method of Lowry.' Bovine serum albumin waa used as standard. 
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